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ABSTRACT 


The problem of allocating different types of jobs among several 
classes of military manpower is becoming increasingly more important as 
military jobs grow more specialized. A model was proposed by Richard C. 
Grinold which constructs a personnel] inventory by rank for each of 
several classes of manpower and then allocates that inventory to meet 
billet requirements. The model is designed for long-range planning 
purposes and produces possible inventory based on an optimization scheme 
that sets permitted errors in the allocation. The thesis presents a 
review and demonstration of the model based on the U.S. Navy officer 
corps, a discussion of implementation considerations, and further work on 


optimization schemes. 
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I. INTRODUCTION 


The problem of allocating different types of jobs among several 
classes of military manpower is becoming increasingly more important as 
military jobs become more specialized. A model was proposed by Grinold 
[Ref. 1] which constructs a personnel inventory by rank for each of 
several classes of manpower and then allocates that inventory to meet 
billet requirements. The model is designed for long-range planning 
purposes and produces possible inventories based on an optimization 
scheme that sets permitted errors in the allocations. 

The next section reviews the model and offers several points to 
consider concerning its implementation based on practical problems during 
the preparation of this thesis. The review highlights the mathematical 
foundations of the model and defines the essential variables and 
parameters. A complete discussion of the mathematical concepts is found 
in Grinold [Ref. 1]. 

Section III presents a means for measuring the error between require- 
ments and inventories and then expresses this measure as the objective of 
3 types of mathematical programs: linear (goal) programming, elastic 
Programming, and quadratic programming. Sections IV, V, and VI present 
a variety of possible formulations for the three programming models from 
which analysis may be accomplished. The formulations presented include 
those proposed by Grinold [Ref. 1] plus several new ones. For the 
quadratic programming models, analytical solutions are derived from which 


computer programs in the APL language were written. 
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Section VII is a demonstration of the model using both APL programs 
and a linear programming optimizer available at the Naval Postgraduate 
School. The demonstration is based on typical data from the U.S. Navy 
officer corps. The last section summarizes the results of the investigation 


and proposes areas for further studies. 


et 





Pl] SREVIEW OF THE MODEL 


Bee PURPOSE 

The purpose of the model is to allocate different types of jobs 
among several classes of manpower in a steady state system. The term 
"allocation" refers to a method for distributing jobs among several 
classes of people or people among different types of jobs. The 
allocations are determined by arrays of sharing fractions which are 
defined below. These sharing fractions represent an actual or possible 
management policy. 

The term “steady state" refers to the equilibrium condition prevailing 
in the model. This condition arises from the assumption of a fixed 
longitudinal manpower flow in each manpower class used for the analysis. 
This concept is defined for the model in Grinold [Ref. 1], and is discussed 
in detail in Grinold and Marshall [Ref. 2], and in Bartholomew and Forbes 
(Ref. 3]. Further discussion of the steady state conditions is found in 
the section which defines the steady state variables. 

The final item of interest in stating the purpose deals with the 
planning aspects of the model. Since this analysis is done for steady 
State conditions, it can be used most effectively to test current or 
alternate management policies for their lona-term effect and provide 
answers to the "what if" policy questions. As an aggregate planning 
tool, the model is intended to test long-term policies that may be used 
later for daily operation of a manpower system. Since conditions of 


Steady state may be significantly different from current conditions, 
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the analyst should be aware of the limitations of the model. Often, 

the relative differences between competing policies or the magnitude of 
the effect of a policy may be the point of interest for which an analysis 
is being performed. The model should be used with these considerations 
in mind. 

The model is motivated by a study of the U.S. Navy officer corps 
and examples from this manpower system will be used when explaining 
concepts as an aid to better understanding of the model. A description 
of the structure of the U.S. Navy officer corps for the unrestricted line 


officers is found in NAVPERS 15197A [Ref. 4]. 


B. DEFINITIONS AND TERMINOLOGY 
1. Attributes 
The definition and understanding of variables in the model is 
complicated by the large number of variables present. The most important 
variables relating to concepts are presented in this section. Other 
variables relating to specific optimization techniques or used for 
mathematical clarity will be introduced as the need for them arises. A 
list of variables used in the model is found in Appendix A. 
The requirements and the inventories possess three common attributes. 
They are: 
e@ MANPOWER CLASSES--the types of manpower; indexed by k=1,2,...,K. 
@ CAREER STAGES--the experience levels required; indexed by 
ete creer Ls 
e JOB CLASSES--the types of jobs to be filled; indexed by j=l,2,...,Jd. 
The clear definition of attributes is closely related to the purpose 
of analysis. In the demonstration which follows, the manpower classes 
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represent types of officers (e.g., Pilot, Naval Flight Officer); the 
career stages represent ranks (e.g., Lieutenant, Commander); and the job 
classes are naval officer billet categories (e.g., 1310, 1320). 

The level of detail is also closely related to the purpose 
of analysis. For example, the manpower class of pilot could be further 
broken down into pilots by type of aircraft, but this imposes increased 
data requirements which may not be needed if the purpose of analysis is 
to study Navy-wide trends rather than issues specific to the Navy's 
aviation community. 

The attributes may apply to either requirements or inventories. 
In Navy terminology, requirements are usually stated in terms of "billets." 
Some confusion may arise since this is synonymous with the term "job." 
The reader should bear in mind that requirements may be referred to as 
billets or jobs which are distinct from the attribute of "job type." 
Likewise, the term "inventories" is used interchangeably with stocks, 
personnel, and people, but is distinctly different from the term "manpower 
types." The use of these terms should be clear from the context in which 
mney OCCur. 

The demonstration highlights another problem found in practical 
application. The types of people and types of jobs are often expressed 
by the same terminology. For example, 1310 is both an officer designator 
and a billet designator for the manpower class "pilot" and for the job 
class "pilot." To avoid confusion, the word equivalent "pilot" will be 
used when speaking of 1310 officers and the billet designator "1310" will 


be used when speaking of a corresponding type of job. 
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Other interpretations of the model based on current military 
personnel problems are possible. For example, manpower classes could be 
defined as enlisted personnel by mental group category and the purpose of 
analysis could be to study the effects on billets of increasing require- 
ments for the higher mental groups with a shrinking manpower pool from 
which to provide accessions. Another possibility would be to define 
manpower classes by sex and job classes by types of sea or shore billets 
and then study the implications of increased accession of women on 
sea-shore rotation cycles. 

2. Steady State Variables 

Let n. be the number of time periods an individual must remain 
in the organization before entering stage (i+l1). For example, naval 
officers in the rank of lieutenant (jg) will remain in the Navy four 
years before being promoted to the rank of lieutenant, thus n. for 
lieutenant (jg) is four. The total amount of time a person may spend in 
stage i is: 

i eee eens 
An example of the amount of time spent in each career stage of the Navy 
officer corps will be given in Section VI (Fig. 4) for the demonstration. 
The steady state variables are defined as follows: 
e Survivor fractions by rank and manpower class--this is the matrix 
S indexed over (I,K). The (et element, Sis represents the 
fraction of type k people who remain in the organization more than 

1 stages. The element Si, may also be interpreted as the 

probability that a type k person remains in the organization 


longer than n. time periods. 


iS 





ye element 


e Expected stage length is the matrix W whose (i,k 
represents the amount of time a person in manpower class k is 
expected to spend in stage i. W may also be defined as the 
product of the probability of a type k person reaching stage i and 
the length of stage i: 

Wig = SaqhMy > 14-4): 
A more complete discussion of survivor fractions and expected stage 
lengths may be found in Grinold [Ref. 1] and Grinold and Marshall [Ref. 2]. 
3. Requirements 
Requirements are expressed in three ways: 

e Billets by rank and job type--the number of jobs at each stage to 
be filled. This is a matrix B whose elements are indexed by 
(1,3). 

e Billets by rank and people type--the number of people required at 
each stage. This is not the normal way in which requirements are 
stated and should not be confused with personnel] inventories. 
This is a matrix P whose elements are indexed by (i,k). 

e Target allocations--the number of (i,k) people desired in (i,j) 
jobs. This 1s a three-dimensional array T with elements indexed 
BK; 5) )- 

4, Inventories and Accessions 

Finding optimal inventories and accessions will normally be the 
goal of this model. They are defined below. 

e Accessions are denoted by a vector Y whose me element represents 


the number of people entering manpower class k each year. 
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Inventories are expressed in three ways: 

e Inventories by rank and job type--the matrix X wnose elements are 
indexed by (i,j) defines the number of jobs at each stage actually 
filled. This is not the normal way in which inventories are 
stated and should not be confused with billet requirements. 

@e Inventories by rank and people type--the actual number of available 
people by manpower class and rank. This is the matrix Z whose 
elements are indexed by (i,k). 

e Actual allocation--the actual number of (i,k) people who will fill 
(i,j) jobs. This is the three-dimensional array A with elements 
indexed by (k,i,j). The notation for requirements and inventories 


is Summarized in Table l. 


eee 
REQUIREMENTS AND INVENTORIES 


Job Type People Type Peoole Type, Rank 


and Rank and Rank and Job Type 
Requirements D5 Pi. EG 
Inventories Xa Zip aj 


5), Sharing Rules 


There are two rules used to distribute the (i,j) jobs among the 
(i,<) manpower types or vice versa. These rules are Known as job sharing 


and people sharing. 


ne 





e Job sharing is defined by a three-dimensional array F whose 
een 
(cece 


be performed by manpower type k. 


element gives the fraction of (i,j) jobs that should 


e People sharing is defined by a three-dimensional array G whose 
ehh 
ene) 


Should perform job type }. 


element gives the fraction of (i,k) manpower types who 


C. DETERMINATION OF REQUIREMENTS 
The target allocation (T), which is the number of type j jobs to be 
filled by class k people of rank i, is determined in one of two ways: 
1. By specifying the desired billets by people type and rank (P) and 
the rule for sharing these people among the different billets (G). 
Note that G times P distributes the requirement for the k~ class 
of people in rank i over the job class j and produces the target 
allocation T: 
Sig Pik ~ "kij’ 
2. By specifying the billets by rank and job type (B) and the rule 
for sharing these billets among the different types of people (F). 
Note that F times B distributes the j type of job in rank 1 


over the people class k and also produces the target allocation 
a 


Kaig Oig = Tkig: 


In order to maintain consistency, the following equation must hold: 


The argument from which this type of equation is derived will be 
referred to as the "consistency rule." Furthermore, all people classes 


must be filling some job and all job classes must be filled by some 


people which means that the equation: 
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Tag = | 
5 Ral 
j J 


implies that 


Pee Nici) ig (2) 
J 
and the equation: 
ioe =k 
kK sad 
implies 
bs 3 > : SK ig Pik (3) 


Thus, one may either start with F and 8 and calculate P from equation (2) 
and then G using equation (1), or start with G and P and calculate 8 from 
equation (3) and then F from equation (1). 

For the Navy example which follows, it should be noted that the usual 
manner of determining target allocations will be to start with a given 
requirement of job types by each rank (8) derived from the force structure. 
Since it is more common to think of people being assigned to jobs rather 
than jobs to people, a typical sharing policy is more easily understood 
in terms of distributing people over the required job types (G). Thus, 
it would be helpful to be able to calculate F and P based on 8 and G. 
Unfortunately, this cannot be done directly. One solution to this 
problem is a least squares approach using equation (3). First determine 
the least squares fit values for P, and then solve for F using equation 


(1). The least squares objective is: 


' \¢ 


ine 8 = ew 
oa ( Ik ig Pik 1a 


foe) Kk 


i 





Notice that this minimization is accomplished by minimizing each of 


the following I expressions: 


Peden > Deve fOr I-1,2,...,;1. 
a | (y kij “ik ae 


For a particular i, the above can be written in matrix notation: 
min (G'P - 8)' (G'P - B). 


Thus, the jth 


column of P is given by the optimal solution: 
> Ae Cle ee 
Now calculate f using equation (1) to obtain: 
al 


kij . 
0 aa D5 = 0 


D. DETERMINATION OF INVENTORIES 

The actual allocation (A), which is the number of people of type k in 
rank i who will be available to fill job type j, is determined using the 
Same consistency rule invoked for equation (1): 


ij Zk * ig 7 Tkig MHZ > (4) 
and the implications similar to equations (2) and (3) above must also 


hold. The equation: 


fg... = | 
oak 
a 
implies that 
Za > eae a) (5) 
J 
and the equation: 
ie | 
k a 
implies 
a5 ; Ikij 7ik (6) 
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In addition, the model 1s constrained by the steady state requirement 
that the expected inventory of type k people in rank i is equal to the 
number of accessions of type k people times the expected duration of type 
k people in rank 1. In equation form, this constraint is: 


aime ak ok (7) 


Since inventories may be computed from accessions (Y) and expected 
stage lengths (W), the optimization techniques to determine inventories 
found in the next sections will usually have Y as the decision variable. 
If the inventories (X or Z) are the decision variables, then a method is 
needed to determine the corresponding accessions. From equation (7), 
note that for each vy there are I equations. Once again, a least 
Squares approach may be used to find the best fit values of Yr for this 
system of equations. The objective is: 


| De 
7 Bae eee ae 


This yields the result: 
-l 


2 
eens ie ctk) (72° 


ik) ONG Kae camera 


on 





ee Ora LON TECHNIQUES 


A. PURPOSE 

This section discusses the use of mathematical programming to deter- 
mine desired values for the variables defined in Section II. The degree 
to which these techniques may be used depends on availability of data, 
level of detail in attribute definition, and desired accuracy of results. 
The scope of the problem to be solved determines the type of the optimi- 
Zation technique and the way in which the problem is formulated. 

The purpose of this section is to present a method by which the 
quality of an allocation may be measured. Such a measure will be used as 
the basis for the objective functions for three types of mathematical 
programs: linear (goal), elastic, and quadratic programming. Formula- 


tions for each of these will be presented in subsequent sections. 


B. FORMULATION OF GOALS 

In formulating the optimization, special care must be taken in 
determining the objective to be achieved and the decision variable to be 
determined. The overall objective is to match actual allocations to 
target allocations with respect to a decision variable. The methods 
described here seek to do this by minimizing the weighted percent differ- 
ence between inventories and requirements. The decision variable to be 
determined may be any one of the previously defined variables. All other 
variables are then treated as parameters. 

Listed below are the possible decision variables and a brief discussion 
of how they may affect the optimization. 


CL 





@ Requirements--Errors may be minimized with either billets by job 
type and rank (B) or billets by people type and rank (P) as the 
decision variable. Normally B is treated as a parameter derived 
from the force structure. For example, if ten aviation squadrons 
of a certain type are required by the force structure, then 
various numbers of officer pilots of different ranks will] be 
needed to fill the aviation type billets found in those squadrons. 
On occasion, it may be of interest to allow the requirements to be 
decision variables and check the resulting optimal values given 
fixed inventories and sharing rules. 

e Inventories and Accessions--The decision variable may be either 
the inventories or accessions. Due to the steady state assumption, 
it is sufficient to know the accessions and the expected length of 
time in each career stage in order to calculate expected inven- 
tories. For this reason, the accessions (Y) will normally be the 
decision variables. This approach requires that billets (B or P), 
Sharing rules (F or G), and expected duration in each stage (W), 
be treated as parameters. The objective functions in Grinold 
[Ref. 1] are of this type. 

@ Expected Duration in Each Career Stage--This variable may be 
the decision variable if billets, accessions and sharing rules are 
Known. If certain assumptions are made about the conditions from 
which an optimal W is derived, then it is possible to relate the 
optimal W to a survivor function. 

@ Sharing Rules--Minimizing errors with F or G as the decision 


variable will result in the optimal "mixing" policy. If the 


ZS 





objective is to find the optimal policy based only on billets, 
inventories and expected stage lengths, then this approach is 


reasonable. 


C. MEASUREMENT OF ERROR 

As stated previously, optimization is performed to minimize the 
weighted error between requirements and inventories. The error for 
the linear models and quadratic models may be expressed as the excess or 


deficit of the three quantities shown in Table 2. 


TABLE 2 
iiPeEs OF ERROR 


Et - ay 44 : the error between targets and 
allocations by job type, rank, 


and people type. 


kKij 


©: 
{ 


eo IZsy - Pei : the error between the inventory 
of (i,k) people and the corres- 
ponding billets. 


e. |b. . - x. «| : the error between the (i,j) billets 
J J J and the corresponding inventory. 


The first error shown involves all three attributes and may be 
obtained from either the second or the third one since these errors may 
be distributed using the job sharing or people sharing rule: 


d = f 


Nig Tik = lteig > Sag! = PRig Pay 
where diy and ej are defined in Table 2. 
For this reason, the last two types of errors, those by rank and job 
type and those by rank and people type, will be used in the following 


discussion. In the linear programming model, it is important to know if 
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the error is a shortfall or overfill. In keeping with the notation in 
Lee [Ref. 5], shortfalls will be denoted by the superscript minus and 
overfills by the superscript plus. 

Weighting factors are used to allow trade-offs between billets of 
differing attributes. The weighting factors are based on total number of 
billets to be filled and the relative importance of billets. A shortfall 
of twenty billets of a particular type will have a different impact 
depending on whether there were 100 or 1,000 billets to be filled. The 
error values will have a common metric if they are converted to percent 
errors. Thus, a better measure of the error would be (z-p)/p or (x-b)/b. 
However, the actual measure may still differ depending on the nature of 
the job type, rank, or personnel type. For example, a five percent 
shortfall of pilots who are of the rank captain may not be as serious as 
a five percent shortfall in pilots who are lieutenants. Thus the percent 
shortfall in meeting a critical target should be weighted more heavily 
than that of a less critical target. [In order to handle this problem, 
additional parameters indicating the percent trade-off between ranks, job 
types, and people types must be introduced to the model. One way in 
which this may be done is to specify parameters which wil] serve as 
benchmarks to measure the permitted range of unit error in overfilling or 
underfilling a requirement. 

For example, let a five percent unit error be the permitted shortfal| 
of 1310 (pilot) lieutenant jobs. Now for another job, say 1310 captains, 
the relative trade-off for an underfill may be found by answering the 
question, "What percent under target in the assignment of 1310 captain 


billets is as serious as a five percent shortfall in filling 1310 
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lieutenant billets?" In this way, essential judgments about the trade-off 
between differing ranks and job types may be made. The same process 
could also be done for ranks and people types or the consistency equation 
could be used to produce the trade-off parameters. The notation for these 
parameters is shown in Table 3. These trade-offs may now be added to the 
penalty function by multiplying the percent error by the inverse of the 
permitted unit error. Thus the actual weighting factors for each error 
is the inverse of the permitted unit percent error times the requirements. 
Let u and v be the weighting factors for job types and people types, 
respectively. As before, let the superscript plus denote weights associ- 
ated with overfills and the superscript minus denote weights associated 
with shortfalls. The notation for these parameters is in Table 4. The 
total penalty value is now defined as the weighting factor times the 
error. TJhis penalty value 1s a measure of error that takes into account 
both the varying sizes of the different requirements and their relative 
importance. The notation for the penalty measure derived thus far is in 
able 5. 

Obtaining valid weights is of particular importance to the model 
Since the final result of the optimization could be very sensitive to 
these parameters. There are several ways in which they can be estimated. 
Expert judges could arbitrarily select either the parameters to be used 
or the range of minimal values for excess and deficit. For example, the 
requirement for a particular job type and rank may be 100. If a permitted 
range of error is 95 to 110, then » = 5 percent and 9 = 10 percent. The 
range of error is a useful concept and the range bounds bb, p and p are 
defined in Table 6. If the lower unit errors are greater than 100 
percent, then it would be possible to have a negative value for the lower 
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TABLE 3 
PERMITTED UNIT PERCENT ERROR 


by job type by people type 
and rank and rank 
Shortfall] ot Sey 
overfil] J O. 
TABLE 4 


WEIGHTING FACTORS 


by job type by people type 
and rank and rank 
- | -] - -l 
; an) = a al 
overfill uss = rr aE See (>. Psy) 
TABLE 5 
WEIGHTED PENALTY MEASURE 
by job type by people type 
and rank and rank 
eee ge=rek =. Qe = ee 
Shortfall eeeaew ee 7 ae eK 
Dit Vij D5 lie ca kK Say Dey 
rae X., - D., + + eek 
overfill Us. @.. 2 Vix diy Ss ee atin Es 
wes is Pik Pik 
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TABLE 6 
ERROR RANGE OF REQUIREMENTS 


by job type by people type 
and rank and rank 
lower error bound (1 - ae oar = bi. (1 - 6) es Do 
upper error bound (i> aa b = ar ie d.1) Pap = Day 


bound. To avoid this possibility, it will be required by definition that 
the lower bounds be greater than or equal to zero. It is also possible to 
estimate the permitted unit error by job type if given people type and 
vice versa using the consistency rule resulting in the equations: 

Diy Sig 7 5G, TK MON al eK ans) ¢ 

Diy Sig 7 DS faa A Olead Wl Keel. J 
D. OBJECTIVE FUNCTIONS 

1. General 
Let Ei be the penalty value function of errors in job class j 

and rank i, and let Diy be the penalty value function of errors in 


people class k and rank i. Then three possible mathematical programming 


objective functions would be: 


feels min = ¢f Eig 1.e@., minimizing error by job types; 
I) 

m2) min ££ Divs i.e., minimizing error by people types; 
ie 


fro omin ssf E.. + c Devs i.e., minimizing error by job types 
] 


aaa 
and people types. 


p 
K 


Minimizing the third objective function should produce a close match between 
targets and actual allocations and was used in the mathematical programs in 
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Grinold |Ref. 1]. The use of objective (3) requires more data input 

for parameters than either objectives (1) or (2). Objective functions 
based on each of the three representations will be used in the following 
sections. 

In order to shorten mathematical notation in the following sections 
on specific objective functions, it will be assumed that the decision 
variable of interest is the inventory X and the objective function is 
type (1) from above. Also, the subscripts will be dropped for clarity. 

2. Linear (Goal) Programming 

One objective function that measures the weighted percent errors 
is a piecewise linear convex function that has a value of zero if the 
actual allocation meets the target, and increases at a rate inversely 
proportional to the permitted error in meeting that target. This means 
the function is unity when the decision variable assumes either the lower 
or the upper error bound. Therefore, the function will satisfy the 


following three conditions: 


et ) Becies a texs= bo, 
it ) ee) eX 00 rg 
Ghae Boe ite 


Suem a function, as shown in Figure 1, is: 


xD), (b - x) 


E(x) 3b oe 


Max 


= max (ue, ue’) 
Notice that the penalty value will be relatively small (less than 1) for 
values between b and b and that it will increase at a rate l/(b ¥) if 


fom or 1/(b 8) if b < x. 


co 





| ny 


ee i eee ee 


b 5 B 


Fig. 1. Piecewise Function for Linear Programming. 


3. Elastic Programming 


For this type of programming model, the definition of the error 
will be modified slightly as shown in Table 7. The weighting factors 
will be changed since the actual error is now slightly different due to 
the use of the error bounds rather than the actual requirements. In 
order to notationally indicate this difference, the use of underscore and 
overscore for weights and errors will be used in place of superscript 
plus and minus. This type of programming can use a piecewise linear 


objective function of the following type: 


(b - x) ifx<b. 
vb - 

E(x) = 0 ice ele 
(x_- 5) ibe <x 


eb 
etiam ene Oil) €) 
This function is shown in Figure 2. Notice that the penalty is 0 for 
values of x within the range of permitted errors (b, b) and increases at 
arate 1/(b ¥) if x is less than b or 1/(b ®) if x is greater than 5. Also, 
notice that if b = b = 5 then this objective is the same as the goal 
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TABLE 7 
ERROR FOR ELASTIC PROGRAMMING 


by job type by people type 

and rank and rank 
shortfall Sie ese ee oie Fak 
overfill ej = Xij - D4 diy Peay 


rr ee ee ee eee eee eee eee eee 


b b 


Fig. 2. Piecewise Function for Elastic Programming 


programming objective. Thus, the goal program is a special case of the 
elastic proaram. 

4. Quadratic Programming 

Several quadratic programming formulations are presented in 

section VI. These formulations are easy to use and give quick results 
for formulations that can be expressed without inequality constraints. 
The quadratic objective function to be used must have the same properties 
mentioned for the piecewise linear function and it must also be areater 


than zero. these properties may be stated mathematically as: 


eu 





(7) Begone it x= <p), 
(ean) et xX = bo, 
ein) E(x) = lif x=b. 
(iv) Ey) 2 0; 


Unfortunately, quadratic functions which use the previously defined 
penalty measure must have symmetric penalties in order to satisfy the four 
conditions. That is, » must equal @ and 6 must equal » when considering 
the decision variable z. The following objective satisfies these 


conditions: 





E(x) = E ee) when 8 =. 


This type of penalty function is shown in Figure 3. An approximation 
that is useful if 68 and wy are close but not equal is described in 


Grinold [Ref. 1]. 


| 
| | 
| | 
| | 
b b D 


Fig. 3. Function for Quadratic Programming. 


E. PROGRAMMING NOMENCLATURE 
[It will be useful to have a notational scheme to identify the 
Various programming types, decision variables, and objective functions 
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found in the following sections. Let LP, EP, and QP refer to linear 
(goal), elastic, and quadratic programming, respectively. Let the 
objective function that sums errors by job class be a type 1 objective 
function, the objective function that sums the errors by people class be 
type 2, and the one that sums errors by both job class and people class 
be type 3. Thus, the mathematical program LP1(Y) is a linear program 
that uses the first objective function with Y as a decision variable. 
Similarly, QP2(Y,X) is a quadratic program with the second type of 


objective function that uses Y and X as the decision variables. 


F. ADDING COST CONSIDERATIONS 

Grinold [Ref. 1] also discusses a way in which cost considerations 
may be added to the model. Let Cy be the annualized cost of a type k 
accession. Then the multi-attribute objective function: 


D 


rir C. Yt (1 - a) fer E:, ik} 


+ = 
oe 1 
eg J 1 


z 
k 
where 0 < A < 1, will minimize the objective with respect to both cost 
and weighted penalties. This objective leads to a family of objective 
functions depending on the value of the parameter ». Notice, however, 
that the units of the objective function are now in terms of cost units 
and percent error units. 

The use of such multi-attribute functions is advantageous when 
the annualized cost of one type of manpower is significantly different 
than another type. For example, the annualized cost of a general unre- 
stricted line officer is vastly less than that of a pilot due to training 


costs, pay bonuses, etc. 
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IV. LINEAR (GOAL) PROGRAMMING FORMULATION 


A. GENERAL 

A linear programming formulation for minimization of errors is 
found in Grinold [Ref. 1] and is derived from the goal programming 
concept as described in Lee [Ref. 5] and Hillier and Lieberman [Ref. 6]. 
Goal programming is a modification and extension of linear programming 
that allows a simultaneous solution of a system of complex objectives 
rather than a single objective. Notice that the goal of minimizing 
errors by differing ranks, job types, and people types constitutes a 
system of conflicting objectives. A solution which forces exact achieve- 
ment of goals for the lower ranks may result in extreme deficits in the 
higher ranks while one which exactly meets requirements in the higher 
ranks may result in large excesses in the lower ranks. Thus, the goal 
approach is used to seek the best possible solution given the stated 
trade-offs. 

In implementation of the following programs, the user should be aware 
of the tactical problems in solving LP's. It may be advantageous in 
some cases to solve the dual formulations of the problems depending on 
the number of variations and the number of constraints. An example of 
forming the dual for LP3(Y) and LP3(X,Y) is in Grinold [Ref. 1] and a 


general discussion of the topic is in Luenberger [Ref. 7]. 


B. ACCESSION AND INVENTORIES 
The first formulation will be designated LP1(Y) and is shown in 
Figure 4. Recall that LPl formulations involve errors of rank and 
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job type. The constraint arises from the conservation of flow derived 
from equation (6) and sharing property in equation (7). An alternate 
formulation that considers both types of errors is LP3(Y) as shown in 


Figure 5. This is labelled in Grinold [Ref. 1] as LP-P. 


min as Dus er +07. e- ,|* ey Lvs Oia a Ve Ge 
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Figure 6 displays the formulation which uses both X and Y as the decision 


variables. It is LP3(X,Y) and is found in Grinold [Ref.-1] as LP-B. 


min yy d|a+ et + up e,|* a Lv dt) + > a= | 
i gjlhij kL ik 
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Sets xe fe. —~ ef = bd, £xfor all i,j 
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oN for all i,k 
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for all k,i,j 
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@emexPcCiED STAGE LENGTHS 
Another application using goal programming would be to minimize 

the errors with respect to expected stage lengths (W). Recall that the 
Stage length for rank i, person type k is the expected amount of time a 
recruit in class k will spend in rank i. An added constraint for this 
formulation is that the stage lengths have upper and lower bounds. In 
the Navy officer example, such bounds would reflect navy policy or legal 
restrictions imposed by Congress. Let these upper and lower bounds be 
defined by matrices W and W, respectively. LP1(W) may be expressed as 


Shown in Figure 7. 
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Fig. 7. LP1(W). 


D. SHARING RULES 

Of particular interest is a mathematical program to specify the 
sharing fraction that represents the policy for distributing job classes 
among the people types (F) or people classes among job types (G). An 
additional constraint on this formulation is that F must sum to one over 
the index k and G must sum to one over the index j. Also, the sharing 
fractions may be bounded above and below. Using the same convention as 
used in the case of stage lengths, let F and G be the lower bounds and 
let F and G be the upper bounds on F and G, respectively. Notice that if 
no sharing is possible, specifying that both the upper and lower bound be 
zero sets the value at zero. LP1(G) and LP2(F) are formulated as shown 


in Figures 8 and 9. 
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Fig. 8. LPI1(G) 
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E. INTERPRETATION FROM LP THEORY 


B and d_ represent logical 


In classical LP terms, the errors e, e, d 
variables introduced to maintain the equality constraints. The errors e 
and d may be thought of as artificial variables and the errors e” and d 
as surplus variables and hence, they measure the distance from a particular 
solution to the requirement. These objective functions are mathematically 
equivalent to a Lagrangian form with cost coefficients of the decision 
variables equal to zero (see Duff [Ref. 8: p. 64]). The weighting factors 
u’, ar ve and v_ represent bounds on the variables of the dual formulation. 
In fact, an interpretation of a particular weighting factor is that it 
defines the penalty cost per unit violation of its associated constraint. 
Therefore, weighting factors which are approximately equal have associated 
constraints which are of the same relative importance, if binding. In 
other words, the incremental penalty for violating each such constraint 


iS approximately the same. Such information should be of use in analyzing 


the results of these programs. 


ee 





V. ELASTIC PROGRAMMING FORMULATION 


A. GENERAL 

The elastic programming model is described in Duff [Ref. 8]. The 
formulation is quite similar to that of the goal program in the last 
section. However, there is now a subtle but important difference. In 
the goal program, the objective was to allocate people as closely as 
possible to an exact number of billets (either B or P). In the elastic 
program, the objective is now to allocate people so as to be within a 
Specified range of the billets to be filled. For this reason, the actual 
error will be indicated by a change in the notation to the use of under- 
score and overscore in place of superscripts plus and minus. 

Before starting the formulations, it is necessary to define a variable 


to represent the slack from the upper bounds (B or P). Let q. and ree 


j 
represent these slack variables for the (i,j) and (i,k) inventories. Also 
recall that u, u, v, v are the weighting factors and e, e, d, d are the 


errors from the desired ranges. 


B. ACCESSION AND INVENTORIES 

The first formulation shown in Figure 10 below is similar to the 
first formulation in the previous chapter and is denoted EP1L(Y). The 
variable Y once again represents an upper bound on the accession vector Y. 
An alternate formulation that considers both type of errors is EP3(Y) 


presented in Figure ll. 
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Fig. 11. €P3(Y) 
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C. EXPECTED STAGE LENGTHS AND SHARING RULES 

Elastic programming may also be used to determine expected stage 
lengths or sharing rules. Only one example, EP2(F) in Figure 12, will be 
presented since the formulations are quite similar to those in the last 
section. Recall that an additional constraint exists for the sharing 
rules since they must sum to one over the appropriate index. An example 


of this is EP2(F) shown below. 
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Fig. 12. €P2(F) 


D. INTERPRETATION FROM LP THEORY 

As was the case for the linear (goal) program, the elastic program 
is mathematically equivalent to a Lagrangian form. The variables e and d 
are artificial and measure the distance from a particular solution to the 
lower error bound. The surplus variables e and d measure the distance 


above the upper error bound and the slack variables q and r measure the 
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distance from a solution to the associated upper bound. The weighting 
factors continue to represent bounds on the dual formulation. 

This formulation has an advantage in that the solution space for 
zero penalty is greatly expanded even if the error range is small. This 
has improtant consequences for the implementation and achievement of 


acceptable solutions for problems with many constraints and variables. 
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VI. QUADRATIC PROGRAMMING FORMULATIONS 


A. GENERAL 

The quadratic programming model described in Grinold [Ref. 1] pro- 
vides an alternate method for optimization of allocations. A general 
discussion of the solution techniques for quadratic programs is found in 
Simmons [Ref. 9]. The following programs are expressed with only equality 
constraints and will yield solutions that are found by analyticaliy 
solving systems of equations. This provides a useful and consistent 
means of determining allocations. However, the use of QP as presented 
here must be weighed against two disadvantages. 

@e The penalty weights must be symmetric for both deficit and excess 

types of errors. 

@ The decision variables are unbounded and, thus, the QP may vield 

impractical solutions such as negative accessions. 

Since much of the data input for the model is subjective or highly 
variable, and the purpose of the model is to explore questions concerning 
policy planning, the magnitudes and relative relationships of the alloca- 
tion are more important than an exact choice of allocation. Thus, the 
use of symmetric penalties may not present significant problems. If the 
parameters and data inputs to the model are fairly consistent with real 
world conditions, then the use of unbounded variables is unlikely to 
produce inconsistent solutions which would adversely affect the result of 
the analysis. The appearance of negative values in the solution is 


likely to indicate misspecification of parameters or the attempt to 
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achieve unrealistic goals. For the following sections, it will be 
assumed that py = 0 and 6= 6. 
The QP is similar to a least squares approach in that the purpose is 
to minimize the sum of squared weighted errors. Appendix B contains a 
brief summary of the resulting quadratic forms used in this section. 
Reindexing the variables makes it possible to conveniently express the QP 
as a system of equations using matrix arithmetic notation. For mathemat- 
ical convenience, the following change of variables may be necessary: 
eee eG (8) 
where m=1,2,...M and M=IJ; 
n=K (i - 1) +k (9) 
where n=l1,2,...,N and N=IK. 
In this way, matrices of size (I,K) may be reindexed to become vectors 
of size (N) and matrices of size (I,J) become vectors of size M. Three 
dimensional arrays may also be reduced to two dimensions. For example, a 
three dimensional array of size (I,J,K) indexed over (i,j,k) may be 
reindexed as a two dimensional array of size (M,K) using equation (8). 
In using matrix notation, it will be necessary on occasion to reform a 
vector into a diagonal matrix. This means that the eo element of the 


th 


vector will occupy the (k,k)~°~ position of the matrix and all off 


diagonal elements will be zero. 


B. ACCESSION AND INVENTORIES 

The first formulation to be considered is QP1(Y). The advantage of 
this formulation is that it requires less data input since the job 
Sharing fractions (F) and requirements by rank and manpower class (P) are 
not needed. Also, once the optimal accession and inventories are 
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determined, values for F and P may be computed based on the optimization 


results. The objective function is: 


C 2 


min EE Us FQ ig Wi Ye > bg } 
4 K 
Now let Ns ak = ae Wey and then reindex (i,j) to (m) and (i,j,k) to 


(m,k) as given above in equation (8). This results in: 


g 


; Z 
ineeoml aan D a) 


m i K 


which, in matrix notation, equals: 


mk Yk 7 


ech oe Uy ecHyY = B) 
where U is now an (M,M) diagonal matrix of the penalty weights. Taking 
the derivative with respect to Y and setting the result equal to zero 
yields the following solution: 
Y = (H'U'UH)! HtU'UB 
Once an optimal value for Y has been determined, inventories may be 
calculated using equations (7) and (6) from Section II. 
The formulation QP3(Y) minimizes the objective: 
min £2 re iV 
ae: in} 
eed 
Subject to the sharing constraint 


2 
el ee ae ea 


a ee 
ij 


E 
K 


Xi5 * . cen ie 
Substituting the constraint into the objective, letting Nik - 9k i j Wes 
and then reindexing from (i,j,k) to (m,k) as given above in equation (8), 
yields in matrix notation: 

inet UitHyY = B) + Y'OY = 2Y'R +c 


where Q is a (K,K) diagonal matrix whose diagonal elements are given by 


S macs 
Tkk © Fak Mik § (10) 
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. _ 2 
R is a (K) vector with elements aes : Vey Day Wer 3 


and c is a scalar constant, namely c= 2 f£ ve rs , 


Taking the derivative and setting it aa to zero results in 

BU Uive= HUB + QY - R= 0. 
Finally, the solution is: 

Y = (H'U'UH + Q]7) [H'U'UB + 2]. 
This model can be found in Grinold [Ref. 1] as UQ-P. An advantage to 
the use of this model over other QP3 formulations is that it uses the 
people sharing rule G which may be more available as a data input than 
the job sharing rule F. 

An alternative program which uses the same objective function is 
QP3(X,Y). In this case, minimize the objective subject to the sharing 
constraint: 
= OenOr ed lil. Kk. 


Weve ST 
k “k 
j 


Both X and Y are now the decision variables and the constraint cannot 


as 


be substituted into the objective function directly. In order to express 
the constraint in matrix notation, it is necessary to define two new 
matrices, F* and W*. First, F* is defined by the elements: 


* er ei St 


Pikig > | 
0 otherwise 

where mele erie kine AK. and. j=1,.30,d. 

Now reindex from (i',k,i,j) to (n,m) as given by equations (8) and (9) 

above to form the properly shaped F*. This results in a block matrix in 


h 


which each of the blocks is a diagonal matrix of the jt column of F 


for a particular (i,k). W* is formed in a similar manner. Let: 
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i ; wi, for k =k! 
une 0 otherwise 

where tk = Fe eee and (alo. «5 1 
Next reindex W* from (i,k',k) to (n,k) according to (9) above. 
The result is a diagonal block matrix in which each of the blocks on 
the diagonal is a column of W. Finally, the constraint can be expressed 
in matrix notation as: 

poke oe wey = Oo. 
The objective function in matrix notation is: 

imeem Ue Xe BY + ¥*OY — 2Y'R +. c 
where U, Q, R, and c are the same as defined for QP3(Y). Now form the 
Lagrangian: 

Coe Ut xXe= By rt YOY = 2Y°R +c - A (F*X = WY) 
where A is the vector of Lagrange multipliers. Taking the derivative 


with respect to X, Y, and A yields the following system of equations: 


2U'UX - 2U°UB - F*' = 0, 
2QY - 2R + We =aOe 
F*X - WRY aoe 


This can be written in partitioned matrix form as: 


po -— ie | | 
alee Ome 1/2)Fe* | | X | | U'UB | 
| ae || | | | 
| wan + wn- + -n--0----- — | — | 
| eet ie | | | 
Pech | ¥ | = | R | 
| | | Wa | | | 
| --=- + --- + ---------- — | | ------ | 
| | || | | | 
| Fe | -We | 0 pale io 
| | || | | | 
L Be. : zi 
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where 0 indicates a zero matrix of appropriate size. The optimal Y, X, 


and X 1S given by: 


eee iy Poe | 
_—s>©6 CO PeUUOeen| =(l/2)F** |] Pa ves | 
| | | | | | | | 
eo? | oases eco 7-7 -o = | ee | 
| | | | | | | | 
eee el Oe (1/2)ue | | OR | 
| | | | | | | | 
ess S | Mops eeesmes a2 >> S>- 25° | Pees eo oe | 
| | | | | | | | 
A ees oN | 0 | eo | 
| | | | | | | | 
= = ie 5 = a 


The program QP3(X,Y) is found in Grinold [Ref. 1] with the name UQ-B. As 
can be seen from the result, QP3(X,Y) has amore complicated solution 
than QP3(Y) and has an additional disadvantage in that the results may 
Slightly violate the consistency rule. Thus, a determination must be 
made whether to use both X and Y from the results or to construct the 
inventories solely on the basis of Y. It has been presented here as an 
illustration of the use of Lagrange multipliers and as a means of checking 
the results from QP3(Y). 

Another means of checking QP3(Y) which insures that the consistency 
rule is met would be to solve QP3(X), use equation (5) to determine z, 
and then use the least squares approach presented in Section II to find 
values for Y. Those values of Y may then be used to reconstruct the 


inventories and insure consistency. 


C. EXPECTED STAGE LENGTHS 
The next formulation that will be discussed is QP2(W). The objective 
fimecion is: 
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Now reindexing from (i,k) to (n) and changing to matrix notation yields: 
min (QW - R)' (QW - R) 
where Q is reformed to be an (n,n) diagonal matrix. The solution to 


mis orogram 1s: 


W = [Q'Q]? gir 


D. ADDING ADDITIONAL EQUALITY CONSTRAINTS 

The technique of Lagrange multipliers presented with the QP3(X,Y) 
model may also be used with additional equality constraints. For example, 
Suppose the total number of accessions to be brought into the system js 
Known but the exact number in each manpower class is not Known. This 


information may be added as an additional equality constraint: 


where r 1S now the total number of recruits to be brought into the 

System each time period. In matrix notation, this constraint becomes: 
ei or = 0 

where « is a vector of k ones. Using the QP1(Y) model as an example, 


form the Lagrangian with this additional constraint to get 


Tee TUUUN(HYs=9B) = » (e'Y = 7) 
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Take the derivatives, set equal to zero, and solve as before to get 
the solution: 
| r ine = 
| | | | 
ie | | 
| | | 
| 
| 
| 
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The advantage of adding an additional constraint is that it can be 
used to force the quadratic mode! to search for solutions that will be 
more related to "real world" problems. This helps to mitigate the 
disadvantage regarding unbounded decision variables. In fact, a range of 
possible values could be specified and then the optimal solution within 
the range could be found using a simple numerical search procedure. 

Since the solution to each problem during a search iteration is explicit 
and the resolution needed is to the nearest whole integer, such a proce- 


dure is easy to implement and should arrive at solutions fairly quickly. 





VII. DEMONSTRATION OF THE MODEL 


A, IMPLEMENTATION 

The model was implemented in its quadratic form on an IBM 3033 
computer at the Naval Postgraduate School using the APL programming 
language as shown in Appendix C. For the linear and elastic programs, 
preprocessors were written in APL to arrange the data in standard MPS 
programming format for implementation on the XS-system optimization 
package. MPS is an international standard and such format is easily 
transferrable to other types of computers. A sample output is included 
in the program listing section. In addition, listings of some of the APL 
program functions used to produce the results in this section are also 
included. The ease with which data arrays can be manipulated in APL 
resulted in concise and flexible program functions which performed the 
optimizations and the input-output formatting of the data. Also written 
were two “user friendly” functions for data input and data display, and a 
master driver function that allows the user to select from a menu of 


options. Appendix C contains a brief summary of these programs. 


B. SYSTEM DESCRIPTION 
1. Attributes 
This demonstration is based on an example from the U.S. Navy 
Officer Corps. The attributes chosen for manpower classes, career stages, 
and job classes are the same as those found in Grinold [Ref. 1] and are 


Shown in Figure 13. 
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MANPOWER CLASSES 


Type Officer Type Description 
1 110x General Unrestricted Line (GURL) 
2 111x/116x Surface Warfare 
g Bx Gey x Submarine Warfare 
4 131x/139x Pilots 
5 132%) 137 x Naval Flight Officer (NFO) 


CAREER STAGE 


Stage Rank Years of Service 
il ENS OR =F 7 
Le LT JG 2- 4 
3 eg 4- 9 
4 LCDR 9 - 14 
2 COR 14 - 19 
6 CAPT 19 - 26 


JOB CLASSES 


Type Billet Designator Description 
l 1000 General, Nonwarfare 
2 1050 General, Warfare 
3 1110/1160 Surface Warfare 
4 207 0 Subsurface Warfare 
5 1310/1390 Palos 
6 1320/1370 Naval Flight Officer 
7 1300 General Aviation 


Fig. 13. Attributes for Example. 


2. Expected Time in Each Rank 
To compute W, it is necessary to have stage lengths and survivor 
fractions. Navy usage of survivor fractions is usually found in the form 
of continuation or retention rates. Current data was obtained from 
Deputy Chief of Naval Operations (MPT), OP-01, with continuation rates 
hor Navy Unrestricted Line. The continuation rates were used directly to 
compute expected length of time spent in each rank. This data is shown in 


Appendix 0. The APL function WAITS uses continuation rates and stage 
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lengths as inputs and produces the matrix W. 


the program GWAITS will 


also accomplish this but uses survivor fractions instead of retention 


rates. The resulting steady state matrix W is shown in Table 8. 


GURL 
ENS Nass, 
LTJG 1.674 
Bl Ze003 
LCDR POG 
eur Cno2e 


CAPT 0.596 


3. Billet Requirements and Sharing Rules 


TABLE 8 


EREChE@eS AGE LENGINS 


SURF 


1.942 
1.700 
2507S 
ls 
0.896 
Ose 


SUB 


1.891 
oc! 
2.040 
1.040 
Oa, 
Olena 


EEO 


2.000 
eS 
Ze 
LZ 
Os 
0.445 


rrr WR Pr 


NFO 


202 
.887 
24] 
329 
.869 
AEE 


Information concerning the requirements by rank and job type was 


also obtained from OP-01 and the array B was constructed. 


(See Table 


9). The next parameter to be specified is the billet sharing array. This 


information was taken from Grinold [Ref. 1] with some slight modifications 


and is displayed in Table 10. 


TABLE 9 


BILLET REQUIREMENTS BY RANK AND JOB TYPE 


1000 
ENS 330 
LTJG 768 
LT 1780 
LCDR 1456 
CDR 1001 
CAPT 571 


1050 


0 
0 
o7S 
a 
468 
462 


EO 


2258 
1706 
2145 
1510 
990 
354 


20 


660 
700 
884 
877 
518 
taZ 


54 


1310 


eS 
1970 
e730 
1818 
710 
0 


20 


8 
es 
1408 

586 

54 
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Once B and F are available, it is possible 
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GURL 


ENS 
LTJG 
LT 


LCDR. 


CDR 
CAPT 


SURF 


ENS 
LTJG 
hh 
LCDR 
CDR 
ear | 


SUB 


ENS 
LTJG 
Ei 
LCDR 
CDR 
CAPT 


PILOT 


ENS 
LTJG 


LCDR 
CDR 
CAPT 


NFO 


ENS 
LTJG 
el 
LCDR 
CDR 
CAPT 


1000 


O73 
0.70 
e70 
a0 
Oy) 
Q.70 


1000 


O.18 
0.15 
Oe 
Om 
Os 
On 


1000 


o799 
0.06 
0.06 
0.06 
0.06 
0.06 


1000 


0.00 
0.50 
050 
0 
0.50 
050 


1000 


0.00 
0.40 
0.40 
0.40 
0.40 
0.40 


BILLET SHARING ARRAY 
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Ooo0edcdo 
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1.00 
550 
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Ueo0 
0.50 
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0.00 
0.10 
Q.20 
0.20 
Ole, 
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1050 


0.00 
OF 10 
0.18 
O13 
0.18 
0.18 


1050 


0.00 
OQ. 10 
0.12 
OeiZ 
0.12 
Oe 
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to calculate the corresponding requirements in terms of people types and 
ranks, and the corresponding people-sharing array. The results for the 
requirements by people types is shown in Table 11. An example of the 


people-sharing array is shown in Table 12 for the manpower class of 


pilots. 
TABLE 11 
BILLET REQUIREMENTS BY RANK AND MANPOWER CLASS 
GURL SURF SUB BE LUuOn NFO 
ENS 24] 2317 690 WEB ls 55 
LTJG 538 1821 746 2010 2 
el 1246 2601 1066 4296 1795 
LCDR 1019 2014 1079 2402 O21 
CDR 701 1374 672 1362 541 
CAPT 400 671 279 S52 260 
MASEE, 12 
PEOPLE-SHARING ARRAY FOR PILOTS 
1000 1050 1110 1120 1310 1320 1300 
ENS 0.015 0.000 0 OO 2 30 Coe ou0 
LTJG Colao. 000 0 O 0.980 eee oes 
i OrOz 0.016 0 QO 0.880 OPO. 063 
LCDR 0.030 0.043 0 OOROR 7 5/ Oe 70 
CDR iOS ene Obi 0 Ome Osoc) O 0.380 
CAPT OF Ost 0.236 0 QO Q.000 OrEEenoes 


4. Permitted Unit Errors 
The permitted unit errors were determined based on subjective 
judgment. First the matrix for permitted unit errors by rank and job 


type was estimated. This is shown in Table 13. The entries marked as 1 
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TABLE 13 
PERMITTED UNIT ERRORS BY JOB TYPE AND RANK 
1000 1050 PaO 1120 VSL 1320 1300 


ENS 400 900 900 900 400 i 
LT JG 400 l 900 900 900 400 900 
LT 10 8 8 10 8 10 10 
LCDR 10 8 10 15 8 8 10 
CDR 10 8 9 10 8 10 10 
CAPT 10 10 10 10 l 1 10 


correspond to ranks and job types for which there are no billets. 

Giving a very low percentage will help ensure the optimization does not 
seek to place people in these jobs. The corresponding errors by people 
type and rank were generated Dy the use of the consistency rule. The 
unit errors do not necessarily have to meet this requirement, but it was 
felt that the consistency rule would give an approximation as reasonable 
as further use of subjective judgment. The results were checked to insure 
that the numbers were reasonable, and no changes were made. The matrix 
for permitted unit errors by rank and people type is shown in Table 14. 
These last two matrices represent the permitted error for overfilling the 
billets. The matrices for underfilling the billets were done using a 


Similar method. 


C. ACCESSIONS USING QUADRATIC PROGRAMMING 

The first set of programs to be considered are those labelled QP1(Y), 
QP2(Y), and QP3(X,Y). The results of the optimization with respect to 
the accessions are shown in Table 15. The column marked "Obj. Val." is 
the value of the penalty function that measures both types of errors. 
Notice that QP1(Y) which only minimized with respect to errors by job 
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TABLE 14 


PERMITTED UNIT ERRORS BY RANK AND PEQPLE TYPE 


GURL SURF SUB PILOT NFO 
ENS 280 .0 SoU 924 .2 aco. / 416.5 
LTJG Zay.0 960.5 524 .2 G05. 2 Gol. | 
eT 0 5 ape 1.6 lovags 
LCDR del) 1529 Nee 26 ee 
CDR no Os, Wee 6 ae 
CAPT es, iG i246 ze0 bad 

TABLE 15 


ACCESSION USING QUADRATIC PROGRAMMING 


Obj. 
Vl nau SURF SU mee Oil NFO Total 
OP EY) 455.2 | 704 1404 994 1470 332 | 4464 


QP3(Y) Sg086. |e tz 1416 963 1246 272 | 4209 
QP3(Y,X) 400.8 | 704 1380 550 1388 270 | 4298 


type and rank still did well when compared to the other two. In fact, as 
can be seen from the actual accession values, all three programs had 
remarkably close results. Since QP3(Y) did the best in terms of the value 
of the objective function, the results from this optimization will be 
presented in the remainder of this section. 

While close, the results from all three programs have rather high 
penalty values. Recall that the minimum value that can be achieved 
1S zero, and the actual magnitude of the penalty values shows that our 
actual allocation is far off the mark. Another way to gain an apprecia- 
tion of how close the allocation is to the targets is by examining the 
percent error compared to the permitted errors. Listed below in Table 16 
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TABLE 16 
PERCENT ERROR BY RANK AND JOB TYPE 
1000 1050 1110 1120 ESA0 1320 1300 
ENS 352 0 Lg 54 aS -9 0 


LTJG 90 0 G2 22 22 -58 -12 
a 2 -1 13 8 -18 -51 -32 
LCDR -26 -32 -20 -46 -47 -40 -44 
CDR -14 -17 -8 -35 -29 -6 -19 
CAPT 14 31 24 Zs 0 0 47 


are the percent errors from the QP3(Y) program. This table highlights a 
problem that is probably due to the high rate of attrition of the middle 
level officers. Notice that those ranks are underfilled (as indicated by 
the minus sign) while the junior and senior ranks are experiencing an 
overfill. The most serious problem exists in the rank of Lieutenant 
Commander and in the job classes Submarine Warfare Officer and Pilot. 

The model could be made to actually fill those shortage billets by moving 
the permitted unit error in those ranks and job classes closer to zero. 
However, Since the current parameters are already closer to zero than any 
of the others, such changes would only result in worse overfills at the 
high and low end. A better approach would be to revise the sharing rules 


in an attempt to get a better allocation. 
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VIII. SUMMARY 


A. GENERAL 

The purpose of this thesis was to investigate and enhance the mode] 
proposed by Grinold [Ref. 1] with special attention to the optimization 
techniques and considerations involved with implementation. The model 
looks at requirements and inventories of a manpower system by classifying 
them according to people type, rank, and job type. By the use of sharing 
fractions, which represent management policy, it is possible to express 
billets which are usually defined by job type as requirements for certain 
types of people. The same sharing rules also make it possible to express 
inventories of people by the types of jobs filled rather than by the 
types of people filling those jobs. Further, the assumption of steady 
state makes it possible to calculate all inventories from accessions. 
Thus, billets, inventories, and accessions are linked by a unified and 
rather elegant model structure. 

If relative trade-offs between different billets in the form of 
permitted shortfalls or overfills can be made, then it is possible to 
measure the error between the desired billet structure and the steady 
State inventories. Various optimization techniques may then be used to 
examine the long term result of policy decisions regarding the allocation 
of people to jobs and to test alternative policies. 

This thesis has included several areas of research not presently 
found in Grinold [Ref. 1]. Several considerations were discussed concern- 


ing the implementation of the model. These include: level of detail in 


60 





attribute definitions, problems encountered in implementation, acquisition 
of model parameters, intuitive explanation of weighting factors, and 
advantages of various optimization schemes. New programming formulations 
were presented that used accessions, inventories, sharing rules, or 
expected stage lengths as the decision variables. Elastic programming 
was introduced as a generalization of the goal programming technique. 
Explicit solutions were derived for the quadratic programming mode! and 
the idea of additional equality constraints was proposed. Other enhance- 
ment features included the least squares approach to determine the 
requirements by people type and the job sharing rules given the billet 
structure and the people sharing rules, and using least squares to 
determine accessions from a given inventory. 

In conclusion, this mode! allows the policy planner to examine a wide 
range of policy options using only a few simplifying assumptions and a 
modest amount of input data. Many areas of policy planning may be 
illuminated by simply placing a problem in the framework of the mode! 
Since this requires that the analyst consider trade-offs between various 
people and job types and formalize policy in setting the values for the 
Sharing fractions. This ability to examine the aggregate effect of human 
behavior on an organization's billet structure makes this model a useful 


too! for the policy planner. 


B. AREAS FOR FURTHER STUDY 
There are several! areas related to this model where further research 
Should be done. Suggested areas include the following: 
@ extend the model structure to handle dynamic conditions. If 
the current inventory is Known and the continuation rates of 
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personnel can be estimated, then a year-by-year approach could 

be used to examine various short-range policy problems. This 
would eliminate the limitation of steady state assumptions in the 
current model. 

The optimizations rely heavily on reasonable estimates of the 
permitted percent error in failing to fill all billets. More 
research is needed on how to gather, scale, and interpret data 
that could be used by the model with regard to these permitted 
errors. 

More improvement is needed in adding cost considerations to the 
model. The ability to relate any analysis to "the bottom line" 
should enhance the results of such an analysis. The idea of a 
multi-attribute objective function that considers both errors and 
costs was presented in Section III. Another possibility would be 
to develop a means to econometrically express the trade-off 
parameters. 

The usefulness of the model would be extended if the results could 
be expressed in such a manner that the impact of policy decisions 
could be quickly evaluated. One idea is to use the model to 
create ratios or indices that pinpoint the criticality of various 
allocations. For example, if the "optimal" billet structure based 
on projected inventories were calculated, then the ratio of this 
theoretical structure and the actual structure should provide a 
measure of “how critical" certain billets will be in the future. 
Ratios close to one are desirable, and those far away highlight 


the problem areas. 
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e Since the model is based on steady state assumptions, it is 


desirable to be able to relate the results of the model to these 
assumptions. Assuming certain functions that describe the way in 
which people leave the organization could be used to establish 
such a relationship. The impact of changes in continuation rates 
or survivor fractions on the billet structure needs to be directly 
related to the sharing policies that will be required to reduce 
the allocation errors, and conversely, changes in the billet 
Structure need to be analyzed in light of the behavior of people 
in the organization. 

Further model enhancement is needed in order to introduce the 
effects of stage substitution to the model. In military organiza- 
tions, such stage substitutions exist in the form of "selected but 
not yet promoted." The model has assumed that such stage substitu- 
tion (i.e., personnel of one rank filling a job of a higher rank) 
is negligible or has a net effect of zero. The ability to model 
this would add more credibility to the model and subsequent 


analysis. 
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APPENDIX A 

VARIABLE LIST FOR THE MODEL AND PROGRAMS 
Description 

index of the manpower classes. k=1l,2,...,K. 
index of ranks or career stages. i=l,2,...,I. 
index of the job classes. j=l,2,...,d. 
number of people in rank i required to fill class j jobs. 
number of required people in rank i and manpower class k. 
number of (i,k) people required in type (i,j) jobs. 
inventory of people in rank i who are filling job j. 
inventory of people in rank i and manpower class k. 
inventory of (i,k) people who fill type j jobs. 


fraction of (i,j) jobs that should be performed by 


manpower class k. 


The fraction of type (i,k) manpower who should fill job 
type Jj. 
The expected length of time a person in manpower class k 


will spend in rank i. 


The 


fraction of type k people who will remain in the organ- 


1Zation past rank i. 


The 
and 


Tne 


permitted unit error for underfilling job class j 
Rank. 1. 


permitted unit error for overfilling job class j and 


yank). 


The 


permitted unit error for underfilling rank i and man- 


power class k. 


The 


permitted unit error for overfilling the rank i and 


manpower class k. 


The 


lower bound of permitted error for job class j and 


rank ic 
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etc) 
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The upper bound of permitted error for job class j and 
rank 1. 


The lower bound of permitted error for manpower class k and 
rank 1. 


The upper bound of permitted error for manpower class k and 
rank 1. 


An indicator that is 1 if manpower class k can fill job j; 
O otherwise. 


The error in meeting (i,g)h" requirement. 
The error in meeting (isk requirement. 
The weighting factor for an (i,j) error. 


The weighting factor for an (i,k) error. 
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APPENDIX B 


QUADRATIC FORMS FOR THE MODEL 


1. The General Form 
The quadratic objective functions in this thesis have the same general 
quadratic form which can be expressed in summation notation or matrix 
notation, if properly reindexed. Consider the quadratic form with 
mama iles X(1,1), ACI,1), P(1I,1) and constant r: 
ae 
z 
Pas Xs eae pix, tr. (al) 
j ] 
In matrix form this becomes: 
X'QX - 2X'P +r (a2 ) 


where Q is the diagonal matrix: 


2. The Special Form 
If the above general quadratic form has the following characteristics, 


then a more familiar form can be written for it. 


= y = 


and then (al) becomes: 
2 


C 
ms zy oe 
x. ie a.D-x. ; D. 


oe i 
aes) 7} : 
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ineigatrix motation this is: 
oe OX eae A Db + BG 
(DX Be Dh = BD (a4) 


where D is the diagonal matrix: 


dq 
aT 


This special form has several cases which are of interest. 
Case l 


Let a, = 1. Then (a3) becomes: 


E(x; - b;)° (25) 
j 
Or in matrix form: 
(ee Bye 0k == 8) (a6) 
Case 2 
Let D. = a.c. Then (a3) becomes: 
eee (a7) 
} 
YapaDx =ezX 0 DE = €'D'DC 
(DX - DB)' (DX - DB) 
CGawe DD (X = C)\ . (a8) 
Case 3 
Let F = a,c. and ca ik 2k 
pee FZ = X, this form is like (a7) and has similar results: 
2 2 
ba, cs Fe zy) - c.] (a9) 
1 kK 
Sa matrix notation: 
teen OD (FZ = C) . (al0) 
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Function 


QUAD 


QP3XY 


MAKEPG 


MAKEBF 


MAKEFP 


QP3Y 


INVENTORY 


PEOPRCNT 


APPENDIX C 
APL FUNCTIONS 


Uses Computes Syntax and Description 


All variables QUAD. 
Master Driver Program. This 
function calls all others through 
the use of menu options. 


Bemis bee We \ a ee ae QP3XY. 

a, > Calculates “optimal” accessions 
Y, personnel inventory Z, and 
allocation A, given the data, 
objective, and the billet share 
rule F. 


PG Seis MAKEPG. 
Reconciles the global variables 
B, F, and T with the global 
variables P and G. 


Biss h Gn: MAKEBF . 
Reconciles the global variables 
P, G, and T with the global 
variables B and F. 


B 


~ 
G?> 
ww 

we 
=m] 

w 
—{ 


MAKEFP. 

Reconciles the global variables 
P, F, and T with the global 
variables B and G using least 
Squares approach. 


Ope omiaters Nee Noe SN Reo 

W, 9, 0 Calculates "optimal" accessions 
Y, personnel] inventory Z, billet 
Staffing X, and allocation A, 
given the data, and use of the 
people share rule. 


G, Y, W Lore A INVENTORY. 
Calculates inventory Z, X and 
allocation A, where accessions, 
people share rule, and expected 
Stage lengths are given. 


Bd? 100(D-Z)/P PEOPRCNT. 


Calculates the percentage error 
in inventory. Compare with 4. 
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Function Uses Computes syntax and Description 


BILPRCNT 8B, X 


GWAITS ALPHA, S 


WAITS REMenS 

ROUND 

LABIJ EAgL AG. 
Oe Lies 

LABIK ILAB, KLAB, 
Beas 

LABKIJ ILAB, KLAB, 


SECbewUee cs 


DISPIN Ge ies Baek; 
ILAB, KLAB, 
BAG UR se CS 


DISPOUT As es 
ILABS, KLAB, 
eeaB. Ue CS 


100(B - x)/B 


BILPRCNT. 

Calculates the percentage error 
in meeting billet requirements. 
compare with 6. 


GWAITS. 

Given the survivor fractions 
ALPHA and stage definitions S, 
calculates W the expected waiting 
time in each stage. 


WAITS. 

Given the retention rates RET 
and stage definitions S, cal- 
culates W, the expected waiting 
time in each stage. 


ROUND XXX. 
Takes any array and rounds 
elements to integers. 


LABIJ XXX. 

Takes any array of dimension 
(I,J) and labels rows and col- 
umns according to ILAB an JLAB. 
ILAB must have I rows and JLAB 
must have J rows. ODP specifies 
number of decimal places desired. 
CS specifies minimum column 
spacing. 


LABIK XXX. 
Handles any array of dimension 
(I,K) in the same manner as LABIJ. 


LABKIJ XXX. 

Handles any array of dimension 
(K, I, J) in the same manner as 
LABIJ. 


DISPIN. 

User selects from a menu to 
produce the input variable with 
appropriate labels. 


DISPOUT 

User selects from a menu to 
produce the output variable with 
appropriate labels. Also displays 
error arrays from BILPRONIT and 
PEOPRCNT. 
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APPENDIX € 
COMPUTER PROGRAMS 


TECTSELAY iy 
Teor BAYS 
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